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Abstract. One of the major sources of uncertainty in 
model estimates of the global sea-salt aerosol distribution is 
the emission parameterization. We evaluate a new sea-salt 
aerosol life cycle module coupled to the online multiscale 
chemical transport model NUMB/BSC-CTM. We compare 
Syr global simulations using five state-of-the-art sea-salt 
Open-ocean emission schemes with monthly averaged coarse 
aerosol optical depth (AOD) from selected AERONET sun 
photometers, surface concentration measurements from the 
University of Miami’s Ocean Aerosol Network, and mea- 
surements from two NOAA/PMEL cruises (AEROINDOEX 
and ACE1). Model results are highly sensitive to the intro- 
duction of sea-surface-temperature (SST)-dependent emis- 
sions and to the accounting of spume particles production. 
Emission ranges from 3888Tgyr~! to 8114Tgyr7!, life- 
time varies between 7.3h and 11.3h, and the average col- 
umn mass load is between 5.0 Tg and 7.2 Tg. Coarse AOD 
is reproduced with an overall correlation of around 0.5 and 
with normalized biases ranging from +8.8% to +38.8 %. 
Surface concentration is simulated with normalized biases 
ranging from —9.5 % to +28 % and the overall correlation is 
around 0.5. Our results indicate that SST-dependent emission 
schemes improve the overall model performance in repro- 
ducing surface concentrations. On the other hand, they lead 
to an overestimation of the coarse AOD at tropical latitudes, 
although it may be affected by uncertainties in the compar- 
ison due to the use of all-sky model AOD, the treatment of 


water uptake, deposition and optical properties in the model 
and/or an inaccurate size distribution at emission. 





1 Introduction 


Sea salt is one of the most abundant aerosol species glob- 
ally. It perturbs the radiative fluxes directly by interacting 
with shortwave and longwave radiation, and indirectly by 
acting as cloud condensation nuclei (CCN) and thus alter- 
ing marine cloud brightness and lifetime. It also influences 
heterogeneous chemistry mainly over coastal areas (Lewis 
and Schwartz, 2004) and is co-emitted with organic aerosols 
(Tsigaridis et al., 2013). The major uncertainties in the sea- 
salt life cycle are emission (Textor et al., 2006; de Leeuw 
et al., 2011), water uptake (Textor et al., 2006), and depo- 
sition (Textor et al., 2007). Lewis and Schwartz (2004) es- 
timate the total sea-salt emission to vary from 0.3 Tgyr7! 
to 30Tgyr~! and estimates from models involved in the 
AEROCOM project range from 3 Tg yr~! to 18 Tg yr! for 
year 2000 (Textor et al., 2006). These uncertainties may lead 
to differences of a factor of two or more in the simulated 
monthly averaged concentrations among different models, 
and between simulated and observed concentrations (Textor 
et al., 2006). The lack of comprehensive measurement data 
sets hampers evaluation efforts and the improvement of sea- 
salt models and related parameterizations. For a given region 
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and a given time period, only a few coincident measurements 
of surface concentration, aerosol optical depth (AOD), and 
particle-size distribution are available, and a few emission 
and deposition flux estimates at specific sites and tempo- 
ral intervals can be found in literature. Additional difficul- 
ties arise from biases in satellite retrievals, particularly in the 
most important sea-salt production regions (e.g. Jaeglé et al., 
2011). 

Several approaches are typically used to parameterize the 
sea-salt emission process, from semi-empirical combinations 
of whitecap factorization and concentration measurements 
(Monahan et al., 1986; Smith et al., 1993; Smith and Harri- 
son, 1998; Andreas, 1998; Hoppel et al., 2002; Gong, 2003; 
Petelski et al., 2005; Martensson et al., 2003; Clarke et al., 
2006; Caffrey et al., 2006; Jaeglé et al., 2011; Fan and Toon, 
2011), to empirical methods such as the use of concentration 
vertical profiles from aircraft observations (Reid et al., 2001). 
Parameterizations of sea-salt emission fluxes may account 
for different production mechanisms (e.g. bubble bursting, 
spume cutting), which may depend on different meteorolog- 
ical parameters. The most used parameter is wind speed at 
10m (Uj), but there have also been attempts to include de- 
pendencies on sea surface temperature (SST), wave height, 
increasing/decreasing wind, salinity and other parameters. 
Exhaustive reviews of these efforts and their performance 
can be found in Lewis and Schwartz (2004), O’Dowd and 
de Leeuw (2007), de Leeuw et al. (2011), and Grythe et al. 
(2013). The above-mentioned parameterizations are assumed 
for the open ocean. Production in the surf-zone represents an 
additional open issue (de Leeuw et al., 2000). 

The high hygroscopicity of sea-salt requires water uptake 
schemes based on prescribed growth factors (Chin et al., 
2002) or equations (Gerber, 1985; Ghan et al., 2001) or ex- 
plicit calculations of the condensed aerosol water (Vignati 
et al., 2004). However, their performance is hard to assess 
and it remains an open topic for aerosol modeling (Textor 
et al., 2006). 

In this contribution, we investigate the uncertainties as- 
sociated with sea-salt, open-ocean emission schemes. We 
use a new sea-salt module coupled online to the multiscale 
NMMB/BSC Chemical Transport Model (NMMB/BSC- 
CTM) (Pérez et al., 2011; Haustein et al., 2012; Jorba et al., 
2012), developed at the Barcelona Supercomputing Center in 
collaboration with NOAA/National Centers for Environmen- 
tal Prediction (NCEP) and the NASA Goddard Institute for 
Space Studies. Its meteorological core, the Non-hydrostatic 
Multiscale Model (NMMB) (Janjic, 2005; Janjic and Black, 
2007; Janjic et al., 2011; Janjic and Gall, 2012) allows for a 
bridging of the gap between global, regional, and local scales 
by using consistent dynamics and physics formulations. 

In Sects. 2 and 3 we present the modeling system along 
with details of the sea-salt module development. We imple- 
ment five emission parameterizations following the whitecap 
approach, in order to investigate this major source of un- 
certainty. Particular attention is given to the description of 
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spume particles production and to the dependence of emis- 
sions upon SST. In Sect. 4 we present the observational data 
sets used for the evaluation of the global simulations. Re- 
sults and discussion of 5 yr simulations (2002-2006) are pre- 
sented in Sect. 5. Evaluation is performed against cruise data 
from the NOAA/PMEL Laboratory, sea-salt concentration 
monthly climatologies from the University of Miami Ocean 
Aerosol Network, and monthly averaged measurements from 
distributed AERONET sun photometers. Our results are also 
compared with other recent model studies, such as Jaeglé 
et al. (2011) and Tsigaridis et al. (2013). 


2 Modeling background 


The NMMB/BSC-CTM is a fully online chemical trans- 
port model coupling the atmospheric equations of NUMB 
with the gas-phase and aerosol continuity equations of 
BSC-CTM. At the present stage of development aerosol 
species included in the model are dust and _ sea-salt. 
The implementation and evaluation of other globally rel- 
evant aerosols is underway. Details on the dust aerosol 
module and gas-phase module can be found in Pérez 
et al. (2011) and Jorba et al. (2012), respectively. Ac- 
cording to the features of its unified meteorological core, 
NMMB/BSC-CTM was conceived for short- and medium- 
range forecasting for a wide range of spatial scales as 
well as for climate studies (http://www.bsc.es/earth-sciences/ 
mineral-dust/nmmbbsc-dust-forecast). Due to its fully online 
coupling, several feedback processes among gases, aerosol 
particles and radiation are taken into account by the model. In 
particular, the radiative effect of aerosols is considered, while 
cloud—aerosol interactions are neglected at present. The on- 
line coupling of aerosol optical properties and gas-phase pho- 
tolysis reactions is also under development. 


2.1. The NCEP non-hydrostatic multiscale 
model (NMMB) 


The NMMB is the meteorological core of the modeling sys- 
tem, allowing simulations of scales ranging from global to 
large eddy simulations (LES) in global and regional domains. 
The regional NMMB has been used at NCEP as the re- 
gional North American Mesoscale (NAM) model since Oc- 
tober 2011. The global model is formulated on the latitude- 
longitude grid, by applying conservative polar boundary con- 
ditions and polar filtering, slowing down the tendencies of 
basic dynamic variables (Janjic, 2009; Janjic and Gall, 2012). 
Rotated latitude-longitude grids are employed for regional 
simulations in order to obtain more uniform grid distances. In 
both cases, the horizontal discretization is performed on the 
Arakawa B-grid. In the vertical, the general hybrid sigma- 
pressure coordinate (Simmons and Burridge, 1981) is used 
with the Lorenz staggering. The “isotropic” horizontal fi- 
nite volume differencing technique assures the conservation 
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of a number of dynamical and quadratic quantities (among 
these, energy and enstrophy). More details about the numer- 
ical schemes of the NMMB can be found in Janjic (1977, 
1979, 1984, 2003) and in Janjic et al. (2001, 2011). 

A variety of physical schemes are implemented in the 
model. A list of these parameterizations and their respec- 
tive references were presented in Pérez et al. (2011) and fur- 
ther details can be found in Janjic (1990, 1994, 1996, 2001). 
For our purposes, we shortly recall the parameterizations in- 
volved in the sea-salt aerosol cycle, i.e., surface layer, grid- 
scale cloud microphysics, convective adjustment and pre- 
cipitation, and radiation schemes. Boundary layer, and free 
atmosphere turbulence are parameterized using the Mellor- 
Yamada-Janjic (MYJ) turbulence closure scheme (Mellor 
and Yamada, 1982; Janjic, 2001). In the surface layer the 
Monin—Obukhov similarity theory (Monin and Obukhov, 
1954) is applied (Janjic, 1996) in combination with a vis- 
cous sub-layer parameterization over oceans (Janjic, 1994). 
The wind speed at 10 m (Uj9), which is the key parameter of 
sea-salt production schemes is computed consistently with 
the surface layer parameterization. The friction velocity u* 
is computed as the square root of the surface layer vertical 
momentum transport. 

Grid-scale clouds are parameterized with the scheme of 
Ferrier et al. (2002) including 5 prognostic cloud variables. 
The relevant quantities for the coupling with aerosol pro- 
cesses are the mixing ratios of both liquid and ice cloud water 
and their conversion rates to precipitation. The Betts-Miller- 
Janjic convective adjustment scheme (Betts, 1986; Betts 
and Miller, 1986; Janjic, 1994, 2000) is used for sub-grid- 
scale clouds. Using conservational constraints, the convec- 
tive clouds are represented by reference humidity and tem- 
perature profiles. Both water vapor mixing ratio and temper- 
ature are relaxed toward reference values within a convection 
time step. In the case of deep convection, the reference pro- 
files and the relaxation time are governed by the cloud effi- 
ciency E which depends on convective regime. This is a non- 
dimensional parameter obtained as a combination of entropy 
change, precipitation, and mean cloud temperature (Janjic, 
1994, 2000). The shallow convection parameterization clo- 
sure uses the constraint that the entropy change must be 
non-negative Janjic (1994, 2000). The NMMB uses the op- 
erational Geophysical Fluid Dynamics Laboratory (GFDL) 
radiation package, which includes shortwave (Lacis and 
Hansen, 1974) and longwave (Fels and Schwarzkopf, 1975) 
schemes. Since the coupling with aerosols is not allowed by 
the operational GFDL scheme, the Rapid Radiative Transfer 
Model (RRTM) (Mlawer et al., 1997) was implemented in 
the model (Pérez et al., 2011). By using RRTM, it is possible 
to couple radiation (both long- and shortwave) and aerosols 
by providing aerosol optical depth, asymmetry factor, and 
single-scattering albedo. 
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2.2 The BSC-CTM dust module (BSC-DUST) 


The development of the sea-salt module follows the im- 
plementation of BSC-DUST (Pérez et al., 2011; Haustein 
et al., 2012), i.e., the dust module of NUMB/BSC-CTM. 
BSC-DUST includes 8 transport bins ranging from 0.1 um to 
10 um in dry radius. Within each transport bin a log-normal 
time-invariant sub-bin distribution is assumed. The processes 
considered by the module are dust emission, horizontal and 
vertical advection, horizontal diffusion and vertical transport 
by turbulence and convection, dry deposition and sedimenta- 
tion, and wet removal including in- and below- cloud scav- 
enging from grid- and sub-grid scale clouds. Water uptake 
was not considered. Given the strong uncertainties on the ac- 
tivation properties of dust, solubility is obtained by applying 
an intermediate hypothesis between pure hydrophobic and 
pure hydrophilic aerosol. Both global and regional simula- 
tions of dust optical depth have been exhaustively evaluated 
in Pérez et al. (2011) and Haustein et al. (2012). In this con- 
tribution, the model coarse AOD is calculated from the dust 
and sea-salt components allowing the use of AERONET sta- 
tions affected by dust to be included in the evaluation. 


3 The sea-salt module 


Sea-salt is assumed to be externally mixed with dust and the 
continuity equation is solved for 8 prognostic size-sections: 


ge + (v*V)nge = pen) a a ho a E. (1) 


n 


where gx are the sea-salt dry mass mixing-ratios, v is the 
wind velocity, subscript 4 stands for horizontal operator, 
and F; Saad Re F ae represent sea-salt production, 
sink/mixing, and turbulent diffusion terms, respectively. Ad- 
vection and diffusion are analogous to those of moisture in 
NMMB (Janjic, 2009). The production term is detailed in 
Sect. 3.1 and sink processes are described in Sect. 3.2. 

We assume a dry radius lower cutoff of 0.1 um in the size 
distribution. Upper size cutoff values depend on the pro- 
duction parameterization (a detailed discussion are provided 
in Sect. 3.1) and were fixed to 15 um to comprehensively 
account for all the different formation processes. Size-bins 
are described in Table 1. Simulated sea-salt mass and opti- 
cal depth are strongly influenced by the number of size-bins 
adopted, due to the strong dependence of dry deposition upon 
particle size (Witek et al., 2011). Simulated values tend to 
converge above 15 size-bins, while mass loss takes place oth- 
erwise. We employ 8 size-bins which involves a mass loss of 
5 % (Witek et al., 2011) — a negligible quantity compared to 
emission uncertainties — as a trade-off for doubled computa- 
tional efficiency. A sub-bin log-normal approach is assumed 
to calculate different momenta of particle radius, such as dry 
effective radius a =< ie >/< ra > and volume mean ra- 


dius ran =(< we >/< re >)!/3_ We assume the canonical 
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Table 1. Sea-salt size bins and their characteristic radii. rg, tans 


and ee stand for dry radius, dry volume mean radius, and effective 
radius, respectively. 








bin ra(um) = r¥™(um) —_r$!*(um) 
1 0.10-0.18 0.14 0.14 
2 0.18-0.30 0.24 0.24 
3 0.30-0.60 0.43 0.45 
4  0.60-1.00 0.77 0.79 
5 1.00-1.80 1.32 1.36 
6 1.80-3.00 2.27 2.32 
7 3.00-6.00 3.98 4.13 
8  6.00-15.00 7,39 8.64 





log-normal distribution of Lewis and Schwartz (2004), char- 
acterized by a geometric radius at RH = 80 %, a6 = 0.3 um 
and geometric standard deviation o® = 2.8. 


3.1 Emissions 


Strong uncertainties of up to one order of magnitude affect 
the estimates of sea-salt production fluxes. The most widely 
used technique to parameterize sea-salt emission is the so- 
called whitecap method, by which the flux is factorized as a 
product of sea-surface whitecap fraction and production per 
whitecap unit, both terms being affected by significant un- 
certainties. Parameterizations use wind-speed at 10m (Uj), 
SST, atmospheric stability, sea-surface salinity, and ocean 
waves properties (height, age, relative direction respect to 
wind), for which Lewis and Schwartz (2004) and O’ Dowd 
and de Leeuw (2007) provide useful reviews. In this study, 
we implement five widely used whitecap method schemes for 
open-ocean production (surf-zone production is neglected) 
with details provided in Table 2. Labels G03, M86, SM93, 
MAO3, and J11 stand for schemes provided in Gong (2003), 
Monahan et al. (1986), Smith et al. (1993), Martensson et al. 
(2003), and Jaeglé et al. (2011), respectively. G03, M86, and 
SM93 are derived from observational data sets and only de- 
pend on Ujo; MAO3 is derived from laboratory experiments 
and includes SST effects that are size-dependent. J11 emis- 
sions are formulated by multiplying the G03 scheme by a 
SST-dependent function equal for all particle sizes. The func- 
tion was fitted using the GEOS-CHEM model and observa- 
tions. In our work, we keep the function as it was derived 
by Jaeglé et al. (2011). With the exception of SM93, all the 
implemented schemes apply the same wind speed power law 
(U ia) in the whitecap parameterization. Consequently, we 
do not focus on the model sensitivity to changes in this term. 
MAO3 was derived for a temperature interval ranging from 
271 K to 298 K, which does not strictly cover the annual vari- 
ation of global SST. J11 is formulated for temperatures rang- 
ing from 273.15 K to 303.15 K. 

For our comparison, we choose schemes differing in par- 
ticle size and production mechanism description. Figure | 
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Fig. 1. Sea-salt number emission flux at 10m as a function of par- 
ticle dry radius with the different emission schemes (color lines). 
Ujo, rq, and SST stand for wind speed at 10 m, dry radius, and sea- 
surface temperature, respectively. 


shows that the strongest uncertainties appear for the ultrafine 
particles (rq < 0.1 um), which do not play a relevant role in 
the simulation of mass concentration and optical depth, and 
thus are beyond the scope of this work. 

All considered schemes account for sea-salt formation 
from bubble bursting. Spume production is not described in 
M86 and MAO3, while it is represented in SM93 (Fan and 
Toon, 2011), and its treatment in G03 is unclear (and, as a 
consequence, in J11). This leads to significant differences in 
emission fluxes of large particles (Fig. 1). 

In addition, the above parameterizations were merged to 
obtain more comprehensive schemes, such as the combined 
M86/SM3 and MA03/M86/SM93 (Table 2). Hoppel et al. 
(2002) concluded that M86/SM93 may be considered as the 
best candidate to describe sea-salt emissions in the interval 
0.15 um to 15 um in dry radius. M86/SM93 was then ex- 
tended to ultrafine particles in other studies (Caffrey et al., 
2006; Fan and Toon, 2011). 

In this work, we also combined M86/SM93 and MAO3 to 
account for the the SST effect upon sea-salt production. In 
MA03/M86/SM93, MAO3 is applied within its range of va- 
lidity and replaced by M86/SM93 beyond that range (i.e., for 
large particles with rq > 1.4 um). We find a similar attempt 
in the work of Tsyro et al. (2011), where MAO3 is combined 
with M86 (but not with the spume production of SM93). 

We choose an upper cutoff for the particle size around the 
maximum value allowed by the sea-salt production parame- 
terizations implemented in our module. To perform a consis- 
tent comparison, we consider a range of [0.1—15] um in dry 
radius for all the emission schemes, which implies an exten- 
sion of M86, G03, and J11 schemes beyond their formula- 
tion intervals. Because some schemes work with wet radius 
rgg and others with dry radius rg, we assume rgo = 27g to 
obtain emission of dry particles following the water-uptake 
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Table 2. Sea-salt number emission fluxes implemented in NUMB/BSC-CTM. dFy /dr fluxes in units [m-2 s—! um 1), dFy /dlog(r) fluxes 
in units [m~? s~!]; rgg and rg stand for wet radius at RH = 80 % and dry radius in units [tm], respectively. If r is used, dry or wet radius was 
not specified. Uy9 in ms—!. SST in K units. Formul. range stands for the size-range in the original formulation of each parameterization. 
The assumption rgg = 2rg is used to merge wet and dry radius intervals. All schemes are applied in the range rg €[0.1—15] um. 








Production scheme Reference Mechanism Formul. range 
(emi) E 

ae lGo3 = 1.373 - uz fre OG of 0.057738) 191-607 exp(—C(r30)") Gong (2003) bubbles, spume: unclear rgq €[0.07—20] 

14d 

A=4.7(1 + Orgo) 0-0! 7"80 6 = 30, C = (0.433 — log(rgg))/0.433 
(emi) 2 

hg = 1.373- URSA! rg 1 +0.057r 4:05) - 101-19 exp(—B(r80)°) Monahan et al. (1986) bubbles rgo €[0.8-8] 

B = (0.38 — log(rgo)) /0.65 

aren Pea \2 

dea Ism93 = dop=1,2 Ak (U10) exp [-4 In(22) | Smith et al. (1993); Hoppel et al. (2002) spume rgo €[5-30] 


log(A1) = 0.0676U 19 + 2.43, log(Ay) = 0.9590}. 


ry =2.1,rg =9.2 


— 1.476 


drem 


Taga IMA03 = 3.84- 10-8 - Ujg!! - (aj (ra) SST + Bj (ra)) 


0 = en 4%j,¢2ra)® Bj = Dea Bj,¢ 2ra)® 
rq € (0.01, 0.0725) > j=1 


rq € (0.0725, 0.2095) > j =2 
rq € (0.2095, 1.4) > j =3 





~(emi) ~(emi) 
max (5 Ism93 > ar ro) if Ujg > 9 








>(emi) 
aren | a 
dy IMs6sMo3 = 
arem) y 
J: — IM86 if Ujo <9 
aren) 
are™ —i—I|mao3_—ifrg < 1.4 


—i-— |IMao3M86SM93 = 
aren) . 
—a— |Msesm93 if rg > 1.4 


arem™ 7(emi) 
N. = Sty : 
di WIL = “digg 1G03 + (SST) 


y =0.3 +.0.1(SST — 273.15) — 0.0076(SST — 273.15)? + 0.00021 (SST — 273.15)3 


Martensson et al. (2003) 


Combined M86/SM93 


Combined MA03/M86/SM93 


Jaeglé et al. (2011) 


bubbles (SST dependent) rq €[0.01-1.4] 


SST €[271-298] 


bubbles, spume 


bubbles (SST), spume 


rgq €[0.07-20] 
SST € [273.15 — 303.15] 


bubbles (SST), spume (SST) 





treatment (detailed in Sect. 3.2). Mass emission fluxes F (emi) 
are calculated from number fluxes F’ . ml) as 


dFO™ dn 

(emi) N 3 

F = ——— - — dra. 2 
k ar 3 Palg Old (2) 

bin—k 

The emission mechanism is not explicitly coupled with the 

viscous sub-layer of the NMMB. However, the calculation 

of friction velocity and wind speed at 10m depends on the 


viscous sub-layer scheme in the surface layer. 
3.2 Water uptake 


The sea-salt life cycle is strongly affected by water uptake. 
Hygroscopic growth may increase particles’ radii by a factor 
of 4 or more. Following Chin et al. (2002) we introduced pre- 
scribed RH-dependent growth factors @(RH) =rw/ra, de- 
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rived from the Global Aerosol Data Set of Kéepke et al. 
(1997) and the database of d’ Almeida (1991) (Table 3). ry 
and rq are the wet and the dry particle radius, respectively. 
We assume the same factors for any radius-moment repre- 
sentation, such as effective and volume-mean radii. 

Given ¢(RH), the water-uptake process is fully described 
by extending any dry particle parameter to its respective wet 
value. In particular we obtain wet particle radius and density 
as 


ra > tw=O''a (3) 
Pa > Pw = fara+ (1 — fa) Pwater, (4) 


where water is the density of water and fg is the volume 
fraction of dry aerosol ( fa = @~°). The dry sea-salt density 
is assumed pg = 2160kgm~? for every size-bin. By using 
this simplified approach, all aerosol processes affected by 


Atmos. Chem. Phys., 13, 11735-11755, 2013 


11740 


Table 3. Sea-salt hygroscopic growth factors @ = rw/rgq at different 
ambient relative humidity (RH) values. ry and rg stand for dry and 
wet radius, respectively. 


RH(%) ¢ 
< 50 1.0 
50-70 §=1.6 
70-80 1.8 
80-90 2.0 
90-95 2.4 
95-99 2.9 


> 99 4.8 


hygroscopic growth are easily reformulated by extending the 
parameterizations used in the dust module (dry aerosol) to 
the wet-particle case, i.e., by applying Eqs. (3) and (4). In 
the following we present a short review of the parameteri- 
zations used in the aerosol module of NNUMB/BSC-CTM, 
pointing out the extension to wet particles in the sink and 
mixing terms. A more detailed description of each scheme 
can be found in Pérez et al. (2011). When not otherwise spec- 
ified we refer to r}™ as rq for brevity. 


3.3 Deposition and convective mixing 


Sedimentation is governed by the gravitational settling ve- 
locity vg «,(@), calculated for each size-bin k following the 
Stokes—Cunningham approximation. vg «,(¢) depends on the 
particle size and thus on the water uptake. 

The dry deposition velocity vdep,«(@), acting at the bottom 
layer, is parameterized following Zhang et al. (2001). The 
dependence on ¢ is introduced in the surface resistance cal- 
culation, which accounts for particle size and density (Slinn, 
1982). 

Wet scavenging fluxes are parameterized both for grid- 
scale (stratiform) and sub-grid scale (convective) clouds. In- 
cloud scavenging flux is parameterized using a solubility pa- 
rameter ¢€; that is defined as the fraction of aerosol contained 
in cloud which may eventually precipitate. For sea-salt parti- 
cles, €, is obtained from Zakey et al. (2006). Since the values 
found in Zakey et al. (2006) for dust represent an intermedi- 
ate between pure hydrophobic and pure hydrophilic hypothe- 
Sis, We assume Es ¢ = 2€dy,¢. This calculation of €ss 4 is con- 
sistent with the values used in other state-of-the-art models 
(see for ex. the sensitivity study in Fan and Toon, 2011). Be- 
cause small particles are more probable candidates to act as 
cloud condensation nuclei, €ss ; decreases with size (see Ta- 
ble 4). Grid-scale below cloud scavenging is parameterized 
following Slinn (1984) in which capture efficiencies Ex (@) 
depend on the wet radius and density of the aerosol particles. 

For sub-grid (convective) clouds, the scavenging fluxes 
are coupled with the Betts—Miller—Janjic scheme (BMJ) of 
the NMMB. The convective in-cloud scavenging parame- 
terization employs solubility factors €(k,@) as well. After 
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Table 4. Sea-salt solubility factors és , for each size-bin k. 
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the in-cloud scavenging, the remaining sea-salt is assumed 
vertically mixed by performing a conservative relaxation to- 
wards reference profiles. The parameterization of sub-grid 
below-cloud scavenging is analogous to the case of grid-scale 
clouds. Within shallow non-precipitating convective clouds 
sea-salt is homogeneously mixed within the cloud. 


3.4 AOD calculation 


In order to calculate the sea-salt optical depth, extinction ef- 
ficiencies Oe. are computed with the Mie-theory solving al- 
gorithm of Mishchenko et al. (2002) for each size-bin k and 
each RH range (Fig. 2). Spherical homogeneous particles are 
assumed. The refractive indices were derived from the Global 
Aerosol Data Set (GADS) (K6epke et al., 1997). Extinction 
efficiencies also depend on the sub-bin log-normal geomet- 
ric parameters r® and o® (Lewis and Schwartz, 2004). The 
optical depth is obtained as 


Tak = BaeMak, (5) 


where Max is the layer dry mass loading of each bin and 
Bx,x iS a mass extinction coefficient which accounts for water 
uptake: 


: 307 75 Gn) 
TN ae 
4reit fab) ba,k 


The algorithm of Mishchenko et al. (2002) also provides 
single-scattering albedo and asymmetry factor for radiative 
calculations. The total sea-salt optical depth is equal to the 
sum over all bins: 


(6) 


ar= >) tie (7) 


k=1,8 


The coarse sea-salt optical depth is calculated with a lower 
cutoff value of 0.6 um (the AERONET submicron cutoff) of 
the wet particle radius. In our description, this value is equiv- 
alent to a lower cutoff of the dry particle radius rg given by 


74 = 0.6um/(RH). (8) 
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At maritime atmospheric conditions (RH = 80 %), the sub- 
micron bins significantly contribute to the coarse AOD. Sub- 
bin contributions to the coarse optical depth are calculated 
assuming the log-normal distribution of Lewis and Schwartz 
(2004). Another useful parameter for model evaluation is the 
resulting AOD (total and coarse) from both sea-salt and dust. 
Because of their external mixing, we assume 


Tss+du,A = Tss,A + Tdu,a> (9) 


where the subscripts ss and du respectively refer to sea-salt 
and dust. 


4 Observational data 


Figure 3 displays the location of measurement sites and 
cruise measurement trajectories used in the model evalua- 
tion. Names and coordinates of the sites are listed in Table 5. 
Quantities evaluated are sea-salt surface concentrations and 
AOD. For the station data we use monthly climatologies to 
compare with our simulated 5 yr period (2002-2006). 

We consider AERONET Sun photometer measurements as 
the reference to evaluate the modeled sea-salt AOD. Even if 
algorithms tend to minimize biases due to cloud cover and 
other effects (e.g. Zhang and Reid, 2006), estimates from 
satellites remain highly uncertain and are not used in this 
contribution. Satellite overestimation can reach up to 0.07 
in island stations compared to monthly AERONET-derived 
AOD (Jaeglé et al., 2011). At certain latitudes, the bias be- 
tween satellite and ship AOD measurements may range from 
—0.2 to +0.2 (Smirnov et al., 2011). These biases exceed the 
typical sea-salt AOD value in the remote marine environment 
(~ 0.07, see Smirnov et al., 2011). AOD measurements from 
the AERONET Maritime Aerosol Network (MAN) are not 
used in this work because of complexities in disentangling 
sea-salt and dust contributions from other aerosol species 
(such as carbonaceous and sulfate aerosols) that are currently 
neglected in our model. 


4.1 NOAA/PMEL cruises 


Sea-salt cruise measurements are considered, specifically ion 
concentrations from two cruises of the NOAA Pacific Ma- 
rine Environmental Laboratory (PMEL): the AEROSOLS99 
and INDian Ocean EXperiment (AEROINDOEX) in 1999 
spanning the Atlantic and the Indian Oceans and the first 
Aerosol Characterization Experiment (ACE1) in 1995 cross- 
ing the Pacific Ocean. Concentrations of both Na* and Cl~ 
were measured by ion chromatography (Quinn et al., 1998) 
at 18 m above the sea surface. The experimental aerodynamic 
cutoff diameter was 10 um for all cruises. Instruments were 
kept at constant RH values during measurements. Based on 
these values, Jaeglé et al. (2011) assumed a dry radius cut- 
off of 3 um for AEROINDOEX and ACE]. Hence, we use 
the first 6 dry model bins for the comparison. The ACE] and 
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Fig. 2. Sea-salt extinction efficiency (Q***) at 4 = 500nm, depend- 
ing on water-uptake growth factors @ of Table 3. The values are 
plotted as function of the effective radius of each size-bin (reff, see 
Table 1). 
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Fig. 3. Observational data sets used for the model evaluation: blue 
circles refer to AOD measurements from AERONET, red trian- 
gles to surface concentration measurements from the U-MIAMI 
network; red lines refer to cruise measurement trajectories from 
AEROINDOEX (solid line), and ACE1 (dashed line). 


AEROINDOEX data sets also provide wind speed measure- 
ments at 33 m and 14m above sea surface, respectively. 

The spatial scale of the cruise measurements is around 
600km since they were averaged over temporal windows 
ranging from 2 to 24h; mean ship speeds were around 
24kmh7!. 


4.2, U-MIAMI surface concentrations 


The U-MIAMI network supplied aerosol measurements from 
around 35 stations worldwide between the early 1980s and 
1996 (Savoie and Prospero, 1977). Aerosols were collected 
with high-volume filter samplers and different measurement 
protocols were employed depending on the measurement 
site. We use climatologies from 15 stations (Fig. 3 and 
Table 5). These stations grant good data quality and are 
not affected by surf-zone production (J. Prospero, personal 
communication, 2012). The observed sea-salt mass concen- 
trations (ugm~*) were computed as SS = CI” + 1.47Nat 
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Table 5. List of the observational sites used in this work, classified 
by network. 


























Code Database Station lat lon 
AODS500 nm 
1 AERONET Amsterdam Island 37.81°S TISTPE 
2 AERONET Ascension Island 7.98°S 14.41° W 
3 AERONET Azores 38.53°N 28.63° W 
4 AERONET Bermuda 32.37°N 64.70° W 
5 AERONET — Cape San Juan 18.38°N 65.62° W 
6 AERONET  Ceilap-RG 51.60° S 69.32° W 
7 AERONET Coconut Island 21.43°N = 157.79° W 
8 AERONET _ Crozet Island 46.43°S 51.85° E 
9 AERONET Dunedin 45.86° S 170.51° E 
10 AERONET — Guam Island 13.43°N 144.80° E 
11 AERONET La Parguera 17.97°N 67.04° W 
12 AERONET Midway Island 28.21°N = 177.38° W 
13 AERONET Nauru 0.52°S 166.92° E 
14 AERONET _ Reunion Island 20.88° S 55.48° E 
15 AERONET _ Rottnest Island 32.00°S ~=—-115.50° E 
16 AERONET © Tahiti 17.58°S  149.61° W 
SURFACE CONCENTRATIONS 
a U-MIAMI _ Baring Head 41.28°S 174.87°E 
b U-MIAMI Bermuda 32.27°N 64.87° W 
c U-MIAMI — Cape Grim 40.68°S — 144.68° E 
d U-MIAMI Cape Point 34.35°S 18.48° E 
e U-MIAMI _ Chatam Island 34.92°S  176.50° W 
f U-MIAMI Fanning Island 3.92°N  159.33° W 
g U-MIAMI Invercargill 46.43°S — 168.35° E 
h U-MIAMI _ King George Island 62.18° S 58.30° W 
i U-MIAMI Marion Island 46.92°S 37.75° E 
1 U-MIAMI Miami 25.75° N 80.25° W 
m U-MIAMI __ Midway Island 28.22°N = 177.35° W 
n U-MIAMI Oahu 21.33°N  157.70° W 
oO U-MIAMI Palmer 64.77° S 64.05° W 
p U-MIAMI _ Reunion Island 21.17°S 55.83° E 
q U-MIAMI American Samoa 14.25°S — 170.58° W 





following Quinn and Bates (2005) where both Cl~ and 
Nat measurements were available, and as SS = 3.252Na‘, 
where only Nat concentrations were supplied (J. Prospero, 
personal communication, 2012). Since U-MIAMI measure- 
ments are not constrained by an upper cutoff radius, we per- 
form the comparison with the complete set of model bins. 


4.3 > AERONET AOD 


The AEROsol RObotics NETwork (AERONET) provides 
automatic ground-based observations from sun photome- 
ters in a large number of stations around the globe (Hol- 
ben et al., 1998; Smirnov et al., 2000). The accuracy of 
AERONET sun photometers is 0.01 for AOD (Holben et al., 
1998; Smirnov et al., 2000). We considered a set of 16 sea- 
salt-dominated stations as proposed by Jaeglé et al. (2011) 
(Fig. 3 and Table 5). The three requirements fulfilled by the 
stations are sea-salt contributions to the total AOD greater 
than 50 % as predicted by GEOS-CHEM model, availability 
of Level 2 quality-assured data for all the considered time 
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ranges, and at least 3 yr of data supporting the monthly cli- 
matologies. The evaluation is performed against monthly cli- 
matologies of the AOD at 500 nm. In particular, we focus on 
the AOD coarse fraction, therefore limiting the influence of 
fine aerosol species. 


5 Results and discussion 
5.1 Experimental setup 


We performed global simulations between 2002 and 2006 
and additional simulations covering the temporal windows 
of the cruises. The horizontal resolution used is 1° x 1.4°. 
24 vertical layers are employed and the dynamics time step 
is At = 120s. Meteorological conditions are initialized ev- 
ery 24h using the NCEP final analyses (FNL) at 1° x 1° for 
year > 2000 and the NCEP Global Data Assimilation System 
analysis (GDAS) at 2.5° x 2.5° prior to year 2000. A spinup 
of 1 month for sea-salt is assumed at the beginning of each 
simulated period. The model output is taken every 6h to cal- 
culate monthly averages and every | h when comparing with 
cruise observations. 

The five implemented emission schemes are compared 
with comprehensive data sets of observations dispersed over 
the globe. The dust AOD is indicated with the label DU. 

Feedback processes between aerosols and radiation are not 
considered in any of the simulations. 

At each evaluation site we also compare the simulated 
wind speed with a 30 yr climatology (1981-2010) derived 
from the NCEP/NCAR reanalysis data set (Kalnay et al., 
1996) to evaluate the representativeness of our Syr wind 
speed climatology. 


5.2. Global sea-salt distribution and total budgets 


Fig. 4 displays the global distribution of simulated sea-salt 
production, surface concentration, and AOD at 500nm in 
January and August with M86/SM93. We observe a pro- 
nounced asymmetry in the summer-to-winter variation be- 
tween the two hemispheres and four large regions of maxi- 
mum production. The two largest monthly peaks are found 
in regions with enhanced westerlies, i.e., beyond the horse 
latitudes (lat > 30°N and lat < 30°S). Also two local max- 
ima can be observed in correspondence with the trade winds, 
next to the intertropical convergence zone (around 10° N and 
10°S). While sea-salt production at the southern belt only 
moderately changes with season, the northern belt is affected 
by strong variations during the year with increases in boreal 
winter well above +200% with respect to boreal summer. 
It is well known that these seasonal fluctuations are related 
to the asymmetric variation of the global wind speed pattern, 
driven by the variation of the global atmospheric angular mo- 
mentum (Sandwell and Agreen, 1984). 

Sea-salt production and surface concentration over the Pa- 
cific around 10° N is about half the values found at higher 
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Fig. 4. Seasonal regimes of sea-salt emission (left panels), AOD at 500 nm (middle panels), and surface concentration (right panels) with 
M86/SM93. January and August averages of a 5 yr period (2002-2006) are shown. The label emi refers to emission flux, and sconc refers to 


surface concentration. 


latitudes. Yet, the AOD reaches monthly mean values close to 
the global maximum. Because of the seasonal movement of 
the Intertropical Convergence Zone, the region around 10° N 
in the Pacific is characterized — during winter — by infrequent 
precipitation and low wet scavenging rate increasing particle 
lifetime, in contrast to the strong production belts character- 
ized by wet extratropical cyclone activity. The RH-dependent 
particle size and optical properties in the model also play 
a relevant role in determining the AOD peaks close to the 
intertropical convergence zone. 

Surface concentration and AOD maximum values in the 
Arabian Sea during the boreal summer are due to the strong 
southwestern winds of the monsoon circulation. 

Figure 5 displays maps of annual mean sea-salt emis- 
sion, surface concentration, and AOD with the five emis- 
sion schemes. The two maximum production regions beyond 
the horse latitudes are the most sensitive to the choice of 
the emission scheme. GO3 produces the highest concentra- 
tions with peaks above 35 ug m7? in the southern belt and 
over 25ugm~? in the northern belt. Differences in spume 
production representation are clear when comparing the sim- 
ple M86 with M86/SM93, for which the mean concentration 
is enhanced due to wind episodes exceeding the threshold 
Uio > 9ms~!. The relative importance of the production re- 
gions changes if SST effects are included in the emission 
scheme. The SST dependence in MA03/M86/SM93 pro- 
duces a latitudinal modulation of the emission fluxes and sur- 
face concentration with relative enhancement in the tropics 
and reduction elsewhere. This effect is amplified with J11, 
leading to a change in maximum values of surface concentra- 
tion from the high-latitude belts to the tropics. In particular, 
an absolute maximum value above 35 ug m7? is found over 
the Arabian Sea. 

Sea-salt AOD patterns with M86, M86/SM93, and G03 
are very similar. The southern belt dominates with peaks 
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around ~ 0.1. Peak values around ~ 0.06 are found at high 
latitudes and the tropical Pacific. Relevant differences are ob- 
served with MA03/M86/SM93 for which absolute maximum 
values of 0.1 appear next to the intertropical convergence 
zone. These peaks overestimate the maximum AOD from 
ship measurements gathered by Smirnov et al. (2011) both 
in the remote tropical Pacific (0.07 for total AOD at 500 nm) 
and the Indian Ocean (0.06, east of Madagascar). The use 
of the J11 scheme leads to an AOD pattern and peak values 
very similar to MA03/M86/SM93, with an enhancement of 
the SST latitudinal modulation. 

Table 6 lists the annual model budgets from the different 
emission schemes and other recent studies. To achieve a con- 
sistent comparison, we specify values for five size intervals: 
all bins (ALL), fine bins up to 1 um (F1), fine bins up to 
0.5 um (F2), coarse bins from | um to 4um (C1), and coarse 
bins from 0.5 um to 4 um (C2). 

Total emission, burden, and lifetime are sensitive to the 
parameterization of the emission flux. Lifetime ranges from 
7.3ih with G03, which produces large spume particles inde- 
pendently from wind conditions, to 11.3h with the simple 
M86, which neglects spume production. 

Despite the decrease in total emissions, SST-dependent 
schemes lead to an enhancement of sea-salt lifetime, both in 
the case of MA03/M86/SM93 (with respect to M86/SM93) 
and J11 (with respect to G03). This effect was also observed 
in Jaeglé et al. (2011), where lifetime values are close to ours. 
However both fine (F1 and F2) and coarse (C1 and C2) life- 
times significantly increase with J11 compared to GO3, in 
contrast to Jaeglé et al. (2011) and Tsigaridis et al. (2013). 
This may be related to different treatments of water uptake, 
deposition, and particle size distribution in the models. 

With respect to AEROCOM experiments, the major differ- 
ence is found in the wet deposition fraction, which is around 
0.4 in our model and between 0.2 and 0.3 in AEROCOM A 
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Fig. 5. Annual mean values of sea-salt emission (left panels), AOD at 500 nm (middle panels), and surface concentration (right panels), 
depending on the emission scheme (from top to bottom: M86, M86/SM93, G03, MA03/M86/SM93, and J11). Averages of a 5 yr period 
(2002-2006) are shown. The label emi refers to emission flux and sconc refers to surface concentration. 


and B median models. In particular, J11 produces the most 
compatible value with the AEROCOM inter-model variabil- 
ity. 

Our simulated annual mean column mass load (ranging 
from 5.0 Tg to 7.2 Tg) is only slightly larger than the value 
of Jaeglé et al. (2011) and it is close to the AEROCOM Ex- 
periment A median value and about half of Experiment B 
(12.0 Tg) (Textor et al., 2006). 

Emission is very sensitive to the upper size cutoff value 


and ranges from 3888 Tg yr~! to 8114Tgyr7!. 


5.3 Modeled surface concentrations compared 
with cruise data 


Cruise measurements allow a comparison with model at 
timescales of 2-24 h. Each measurement gathered by the ves- 
sels was averaged on space and time, thus simulated val- 
ues may be affected by errors due to the adopted averag- 
ing technique. We remap the original lat/lon grid at resolu- 
tion Ax, Ay = 1°, 1.4° to a coarser resolution (Ax’ =nAx, 
Ay’ =nAy), matching the characteristic spatial length of the 
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cruise under consideration. The number n is defined as the 
smallest integer satisfying the following conditions: 


Voruise *max(Tops) < nAx (10) 
Voruise *max(Tops) < nAy, (11) 


where Veruise is the vessel mean speed during the cruise and 
Tops is the observation duration, which is not constant. In 
this way, the spatial extent of each measurement is repre- 
sented by a single lower resolution grid cell. We use n = 2 for 
AEROINDOEX and n = 4 for ACE1. Model outputs every 
1h are then averaged over each measurement period. Cruise 
trajectories are displayed in Fig. 3. We recall that the values 
shown in this comparison refer to an upper cutoff of 3 um 
in dry radius, thus we investigate the model’s ability to sim- 
ulate the concentration within the first 6 bins. In this case 
the M86/SM93 scheme is equivalent to M86, since the larger 
particles produced by spume cutting are not taken into ac- 
count due to the observational cutoff. 

Figures 6 and 7 show a good overall correlation for 
AEROINDOEX and a lower correlation for ACEI. Re- 
sults are similar to those obtained in Jaeglé et al. (2011) 
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Table 6. Model sea-salt global budgets and lifetimes compared with other recent studies. The label emi stands for total accumulated 
emission of sea-salt mass (Tg yr—); wetfrac = wetdep/(drydep + wetdep) where drydep accounts for both accumulated dry deposition and 
sedimentation (Tg yr!) and wetdep for accumulated wet deposition (Tg yr); <load> is the annual mean column mass load (Tg yr), 
and lifetime =< load >/(drydep + wetdep) (h). All quantities obtained in our work are averaged over the simulation period 2002-2006 and 
indicated by the label o.w. The labels ALL, F1, C1, F2, and C2 refer to different dry radius intervals (um). 








ALL (0.1-15) FI (0.1-1) Cl d-4) F2 (0.1-0.5) C2 (0.54) 

Study emi <load> life wetfrac emi life emi life emi life emi _ilife 
M86 O.w. 3888 5.0 11.3 0.486 499 27.7) 2033 14.7 75 28.7 2457 16.7 
M86/SM93 O.w. 5440 56 8.9 0.467 499 27.7 2172 148 75 28.7 2595 16.8 
G03 O.w. 8114 67 7.3 0.400 372 24.7 2465 122 57 26.1 2781 13.9 
MA03/M86/SM93 O.w. 5419 65 10.4 0.466 266 32.5 1997 172 48 33.6 2215 19.5 
Jil O.w. 6514 72 9.6 0.368 298 34.9 1979 165 46 37.1 2233 19.0 
Mgo! Tsigaridis et al. (2013) 471 32.6 1916 26.6 
G03! Tsigaridis et al. (2013) 357 «32.6 92327 26.9 
yi! Tsigaridis et al. (2013) 310 31.9 2019 26.6 
G03? Jaeglé et al. (2011) 5200 47 79 0.356 67 = =26.2 2533 11.3 
Jil? Jaeglé et al. (2011) 4600 44 84 0.402 59 24.7 2229 12.0 
AEROCOM A median? _ Textor et al. (2007) 3830 6.5 7.2 0.210 
AEROCOM A mean? Textor et al. (2007) 8200 7.9 12.0 0.210 
AEROCOM A stddev? __Textor et al. (2007) 8200 54 7.1 0.122 
AEROCOM B median? _ Textor et al. (2007) 7740 12.0 144 0.282 
AEROCOM B mean? Textor et al. (2007) 7720 12.7 12.0 0.253 
AEROCOM B stddev? _Textor et al. (2007) 231 3.9 1.7 0.114 





'2 modes: rg € [0.1-1] um and rg € [1-4] pm 23 modes: rg € [0.01-0.5] um, rg € [0.5—4] pm, and rg € [4-10] pm 3models participating in AEROCOM experiments use different parameterizations 
and aerosol size representations. 
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Fig. 6. From top to bottom: simulated sea-salt surface concentration (sconc) with M86 (blue), G03 (green), MA03/M86/SM93 (violet), and 
J11 (cyan) compared to AEROINDOEX cruise measurements (black squares), simulated wind speed (red line) compared to AEROINDOEX 
measurements (black line) (simulated SST is also shown), and simulated precipitation (red line) compared to AEROINDOEX measurements 
(black line) (simulated accumulated wet deposition is also shown). M86/SM93 is not shown since it is equivalent to M86 for rg < 3um. The 
model 1 h-output surface concentrations are plotted with solid lines and averaged over the measurement times (circles). Skill scores shown 
are correlation (7), mean normalized bias (bias), and mean normalized gross error (g.err.). 
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Fig. 7. From top to bottom: simulated sea-salt surface concentration (sconc) with M86 (blue), G03 (green), MA03/M86/SM93 (violet), and 
J11 (cyan) compared to ACE1 cruise measurements (black squares), simulated wind speed (red line) compared to ACE1 measurements 


(black line) (simulated SST is also shown), and simulated precipitation 
accumulated wet deposition is also shown). M86/SM93 is not shown 


(red line) compared to ACEI measurements (black line) (simulated 
since it is equivalent to M86 for rg < 3um. The model | h-output 


surface concentrations are plotted with solid lines and averaged over the measurement times (circles). Skill scores shown are correlation (r), 


mean normalized bias (bias), and mean normalized gross error (g.err.). 


and Witek et al. (2007) with GEOS-CHEM and NAAPS 
models, respectively. AEROINDOEX is simulated with 
a correlation close to 0.6 with M86, G03, and J11. The 
correlation decreases when employing the SST-dependent 
MA03/M86/SM93 (0.49). The mean normalized bias ranges 
from +0.4 % (MA03/M86/SM93) to —36.8 % (G03), show- 
ing an overall tendency of the model towards underestima- 
tion, mainly due to the misrepresentation of the peak around 
day 25. The mean normalized gross error is around 60 % for 
all cases. 

For ACE1, the correlation is low (around 0.36) regardless 
of the emission scheme applied. The overall bias between 
model and observations ranges from —23.7 % (J11) to 26.3 % 
(M86), with J11 being the only scheme with a negative bias. 
The gross error is around 50 %. Because of the larger mea- 
surement timescales of ACE1, part of the errors may be due 
to the inefficiency of the averaging technique. 

Figures 6 and 7 also display the simulated wind speed 
(at 14m and 33 m), obtained by applying the power-law of 
Hsu et al. (1994) to Uj9, observed wind speed, SST values, 
wet deposition and observed hourly precipitation. The mea- 
surements are 30 min-averages and the model output is taken 
every 1h as an instantaneous value for both cruises. Wind 
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speed is simulated with a high correlation for AEROIN- 
DOEX (0.73) and ACEI! (0.81) and biases (—8.5% and 
—3.7 %) and gross errors (24.2 % and 17.2 %) are rather low. 

The plotted SST values highlight the corrections intro- 
duced by the SST-dependent emissions (MA03/M86/SM93 
and J11) on their parent schemes (M86/SM93 and G03, re- 
spectively). For example, the higher SST values taking place 
during the first 30 days of ACE] produce a marked increase 
in concentration. 

Simulated wet deposition peaks are generally in corre- 
spondence with peaks in observed precipitation, although the 
relative intensity may not always be well captured. Spurious 
precipitations and consequent wet deposition overestimation 
affect the concentration peak around day 25 of the AEROIN- 
DOEX campaign. 

5.4 Model surface concentration at U-MIAMI 
measurement sites 


Simulated surface concentration was compared with ob- 
served climatologies at 15 stations of the U-MIAMI net- 
work (Fig. 8). The discussion of the results is complemented 
by Fig. 9, where simulated wind speed is compared with a 
30 yr wind speed climatology (1981-2010) derived from the 
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Fig. 8. Monthly mean simulated (color lines) and observed (black lines) sea-salt surface concentration at U-MIAMI stations. Simulated 
values cover a 5S yr period (2002-2006). U-MIAMI climatologies include interannual standard deviation bars. J11 interannual standard 
deviation is also shown (shaded grey). The label CLIM stands for climatologies. 


NCEP/NCAR Reanalysis (Kalnay et al., 1996). We use this 
comparison in order to check the representativeness of our 
simulation temporal window (5 yr, 2002-2006) with respect 
to a more comprehensive model climatology. 

Overall, simulated sea-salt concentrations are in good 
agreement with observations. Significant model overestima- 
tions of a factor of two or more are found in Invercargill (g), 
and Marion Island (i), where all schemes are above the ob- 
served mean plus one standard deviation. In these sites, over- 
estimation cannot be attributed to an excess of wind speed 
compared with the NCEP/NCAR climatology (Fig. 9). A 
similar overestimation is found in Jaeglé et al. (2011) and 
Tsigaridis et al. (2013) in Marion Island (i) and in Inver- 
cargill (g) in Tsigaridis et al. (2013). Both studies use as well 
global models with an horizontal resolution greater than | 
degree. Since these sites are located in regions characterized 
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by complex topography, we hypothesize that errors may be 
due in part to the low model resolution used. 

In the Antarctic region, schemes show opposite perfor- 
mances in two stations close to each other. In Palmer (0), 
all schemes overestimate surface concentration with the ex- 
ception of J11, which reproduces well the climatology. The 
overestimation of the wind speed with respect to the NCEP 
climatology does not entirely explain the behavior of M86, 
M86/SM93, G03, and MA03/M86/SM93. On the other hand, 
J11 leads to an underestimation of the observed climatology 
at King George Island (h) which cannot be attributed to wind 
speed. Contrasting results are found in this region when com- 
paring our work with Jaeglé et al. (2011) and Tsigaridis et al. 
(2013). 

Significant underestimation is found in Fanning Island 
(f) and American Samoa (q), both located in the tropical 
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Fig. 9. Monthly mean simulated surface wind speed (2002—2006) (blue) and surface wind climatologies (1981-2010) from NCEP/NCAR 
Reanalysis (black) at U-MIAMI stations. Interannual standard deviations are shown for the simulation (shaded grey) and the NCEP/NCAR 


climatology (black bars). 


Pacific. In Fanning Island (f), the low concentration could 
be associated to an underestimation of the wind speed, in 
contrast to American Samoa (q), where the simulated wind 
speed matches the NCEP climatology. Observed concentra- 
tion climatologies present large standard deviations during 
boreal winter in these sites, which is inconsistent with the 
low variability of the simulated or climatological wind speed, 
suggesting a poor representativeness of the observed mean 
concentration. Significant underestimation of the U-MIAMI 
climatological values in Fanning Island (f) and American 
Samoa (q) is also found in Jaeglé et al. (2011) and Tsigaridis 
et al. (2013). 

There is a significant influence of the applied emission 
scheme upon modeled sea-salt surface concentrations. Even 
the introduction of SST-dependence in MA03/M86/SM93, 
which only affects the smaller bins (from 0.1 um up to 1.4 um 
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in dry radius), makes a relevant contribution to the simu- 
lated concentration. Fig. 10 includes scatterplots of observed 
and simulated values (neglecting Invercargill (g) and Mar- 
ion Island (i)) and a scatterplot of simulated wind speed 
and NCEP climatological values. Correlation, normalized 
bias, normalized gross error, and a linear regression fit are 
provided for each scatterplot. G03 generally overestimates 
the climatological monthly mean concentrations, while M86 
underestimates. A significant reduction in bias is obtained 
when using M86/SM93 instead of M86. The best agree- 
ment is obtained with the SST-dependent emission schemes 
MA03/M86/SM93 and J11. Overestimation with G03 may 
be explained by its unclear description of spume particles 
production, as already noted in Fan and Toon (2011). In- 
deed, the emission flux for particles larger than 10 um in dry 
radius is nearly one order of magnitude larger than in the 
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other implemented schemes (Fig. 1). On the other hand, the 
spume production is neglected in M86, which may partly ex- 
plain the underestimated concentration. The introduction of 
spume particles in the combined M86/SM93 improves the 
model results. This improvement is more evident in stations 
and months characterized by frequent episodes of wind speed 
greater than 9ms~!, such as for example during January, 
February and March in Bermuda (b) (Fig. 8). In the tropics, 
where these episodes are infrequent, M86 and M86/SM93 
provide similar results. 

The introduction of the SST-dependence in the emission 
scheme (both for MA03/M86/SM93 and J11) improves the 
overall statistics, with a reduction in bias and gross error 
(Fig. 10). The wind speed scatterplot suggests that biases 
could not be related to a weak representativeness of the sim- 
ulated 5 yr period. 


5.5 Aerosol optical depth 


Simulated coarse AOD are compared with monthly clima- 
tologies at 16 AERONET sites (Fig. 11). The model AOD is 
all-sky in contrast to AOD measurements, which are clear- 
sky. The differences between all-sky and clear-sky results in 
models are currently uncertain and are thought to be moder- 
ate for sea-salt and very low for dust (Shindell et al., 2013). 
To support the analysis, Fig. 12 displays the simulated 
wind speed and the NCEP wind speed climatology at each 
site. Overall, the simulated coarse AOD is in agreement with 
observations. Significant discrepancies are found in Ceilap- 
RG (6), Dunedin (9), Reunion Island (14), and Tahiti (16) 
with all schemes. Overestimation in Ceilap-RG (6) affects 
the entire seasonal cycle mostly due to errors in dust emis- 
sions from South America. Neglecting the dust contribu- 
tion, the nearly constant seasonal cycle and its mean value 
(~ 0.02) are well reproduced. Overestimation in Bermuda 
(4), Dunedin (9), Reunion Island (14), and Tahiti (16) takes 
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place mainly during austral winter (JJA). At Bermuda (4), 
Dunedin (9), and Tahiti (16) model peaks may be partly re- 
lated to an overestimated wind speed (Fig. 12). 

Both at Reunion Island (14) and Dunedin (9) the signifi- 
cant dust contribution leads to uncertainties in the compari- 
son. At Ascension Island (2), Bermuda (4), Cape San Juan 
(5), La Parguera (11), and Midway Island (12), the model’s 
ability to reproduce the dust cycle is decisive for a proper 
simulation of the coarse AOD. 

Results outline a close behavior among G03, M86 
and M86/SM93, in contrast to SST-dependent schemes 
(MA03/M86/SM93 and J11). The latter tend to overesti- 
mate the AOD over warm sea surfaces in/near the tropics 
(e.g. Bermuda (4), Coconut Island (7), Guam Island (10), 
Midway Island (12), Reunion Island (14), and Tahiti (16)). 
These results are mainly related to differences in the emitted 
size-distribution and the hygroscopic growth of sub-micron 
aerosols affecting the coarse AOD. Figure | shows close to 
an order of magnitude difference in the number emission flux 
for particles with dry radius in the range 0.15—1.4um. Hy- 
groscopic growth of particles above 0.6 um adds up to ex- 
plain the higher coarse AOD when using an SST-dependent 
scheme. Simulated wind speeds are stronger than NCEP cli- 
matological winds in Coconut Island (7), Guam Island (10), 
and Tahiti (16) (Fig. 12), which may also partly explain the 
overestimation in these sites. On the contrary, in Bermuda 
(4), Midway Island (12), and Reunion Island (14) the simu- 
lated wind speed is lower than the NCEP climatology. 

In the subset of stations in/near the tropics, Jaeglé 
et al. (2011) found that their model generally underesti- 
mates the observed climatologies and the introduction of 
a SST-dependent term leads to an improvement of their 
model performance. In our work, we also find an increase 
in AOD at low latitudes when applying SST-dependent 
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Fig. 11. Monthly mean simulated sea salt + dust coarse AOD at 500 nm (color lines) and observed coarse AOD at 500 nm (black lines) at se- 
lected AERONET stations. Simulated values refer to averages over the simulated 5 yr period (2002-2006). AERONET climatologies include 
interannual standard deviation bars. J11 interannual standard deviation is also shown (shaded grey). The dust contribution is highlighted with 
the orange line. The label COARSE stands for coarse fraction of AOD; the label CLIM stands for climatologies. 


emission schemes, although it produces an overestimation of 
AERONET climatologies. 

Figure 13 displays scatterplots of simulated and observed 
coarse AOD, excluding Ceilap-RG (6) due to the significant 
errors in the dust component. MA03/M86/SM93 and J11 
produce the largest positive biases (+38.8 % and +27.5 %, 
respectively) and gross errors (above 40 %). The wind speed 
scatterplot indicates that such overestimation cannot be ex- 
plained by a wind speed overestimation. 
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6 Conclusions 


We presented simulations of the sea-salt aerosol global distri- 
bution with the multiscale model NMUMB/BSC-CTM. Since 
the main uncertainties in sea-salt modeling are related to 
the parameterization of emissions, we implemented five dif- 
ferent sea-salt emission schemes and analyzed their perfor- 
mance. We compared global simulations covering the pe- 
riod 2002-2006 with climatologies from the U-MIAMI and 
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Fig. 12. Monthly mean simulated surface wind speed (2002-2006) (blue) and surface wind climatologies (1981-2010) from NCEP/NCAR 
Reanalysis (black) at selected AERONET stations. Interannual standard deviations are shown for the simulation (shaded grey) and the 
NCEP/NCAR climatology (black bars). 


“sea-salt-dominated” stations from the AERONET sun pho- 
tometer network. For the comparison with AERONET, we 


use the coarse fraction of the AOD. 


We found a strong sensitivity of sea-salt aerosol lifetime to 
the emission scheme, ranging from 7.3 h to 11.3 h. The use of 
SST-dependent emission schemes produces an enhancement 


The 


SST 


latitudinal 


modulation 
MA03/M86/SM93 also leads to marked differences in 


and 


the global patterns of surface concentration and AOD 


compared with M86, M86/SM93, and G03. In particular, 
maximum AOD values are reached at high latitudes with 
M86, M86/SM93, and G03, and in the tropics with J11 and 


of the sea-salt lifetime, which increases from 7.3 h with G03 
to 9.6h with J11 and from 8.9h with M86/SM93 to 10.4h 
with MA03/M86/SM93. 
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MA03/M86/SM93. 
SST-dependent emission schemes lead to a clear im- 
provement of the simulated surface concentration, with a 
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Fig. 13. Scatterplots of simulated sea-salt+dust coarse AOD at 500 nm and climatologies from the AERONET network for each emission 
scheme. Ceilap-RG (6) has been excluded, since the dust errors affecting its region. A scatterplot of simulated surface wind speed and 
NCEP/NCAR climatology is provided in the bottom-right panel. The plots are accompanied by y = 21x dashed lines. r, b, and g stand for 
overall mean correlation, normalized bias, and normalized gross error, respectively. A linear regression fit is also shown (y = mx + q). 


significant reduction of bias and gross error. In particular, 
J11 shows the best agreement with observations. However, 
the simulated coarse AOD with J11 and MA03/M86/SM$93 is 
affected by positive biases at several AERONET sites located 
in the tropics. Factors that may explain the AOD overestima- 
tion include the use of all-sky model AOD in the comparison 
and the treatment of the water uptake, deposition, and optical 
properties in the model. Further research may aim at investi- 
gating SST effects upon particle size distribution. 

As in previous studies, the model shows a strong overes- 
timation in sites characterized by steep topography (Inver- 
cargill (g) in New Zealand and Marion Island (i)), indepen- 
dently from the applied emission scheme. Our preliminary 
results with a high-resolution regional simulation suggest 
that smaller scales play a key role in these sites. A detailed 
investigation at these sites is underway. 

The development of the sea-salt module of the 
NMMB/BSC-CTM is a step forward towards an aerosol 
model, including dust (Pérez et al., 2011), black and organic 
carbon, sulfate, and its online coupling with the gas-phase 
chemistry (Jorba et al., 2012) to obtain a unified online 
multiscale chemical weather forecasting system. 
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